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Abstract 

We obtain conditions for the existence and stability of de Sitter attractors in 
the phase space of homogeneous and isotropic cosmology in generalized theories 
of gravity (including non-linear and scalar-tensor theories). These conditions are 
valid for any form of the coupling functions of the theory. Stability with respect 
to inhomogeneous perturbations is analyzed using a covariant and gauge-invariant 
formalism. The relevance for inflationary scenarios of the early universe and for 
quintessence models of the present era is discussed. 



1 Introduction 



In general relativity de Sitter space plays a special role because quantum field theory 
predicts the existence of vacuum energy which is equivalent to a cosmological constant 
A, and the solution of the Einstein field equations with vacuum energy as the only 
material source is de Sitter space. A period of de Sitter-like inflationary expansion of 
the early universe has come to be regarded as the canonical solution to the horizon, 
flatness and monopole problems that plague standard big bang cosmology. As a bonus, 
inflation provides a mechanism for generating density perturbations through quantum 
fluctuations of the inflaton field, seeding the structures observed in the universe today 
[1, 2, 3]. 

In most inflationary models based on general relativity the expansion of the uni- 
verse described by the scale factor a(t) of the Friedmann-Lemaitre-Robertson- Walker 
(hereafter "FLRW" ) metric 

ds 2 = -dt 2 + a 2 {t)(dx 2 + dy 2 + dz 2 ) , (1.1) 

is approximately exponential. This accelerated expansion is achieved if the dynamics of 
the universe are dominated by a scalar field self-interacting through a potential V(<f)) 
that has a plateau such that V(4>) ~ Vb =constant for a certain range of values of 0. 
While evolves through this interval the potential mimics a cosmological constant. The 
corresponding solution of the Einstein equations has the form 

a(t) = a e H{t)t (1.2) 

where 

H(t) = H + H 1 t+ ... (1.3) 

and do, Ho, and Hi are constant, with \Hit\ « H . In other words spacetime is close 
to the de Sitter solution and the scalar field rolls slowly over the plateau of the potential 
("slow-roll approximation" [4, 2]). The dynamics of a scalar field minimally coupled to 
the spacetime curvature are described by the Klein-Gordon equation 

dV 

+ 3F0+ — = 0. (1.4) 

d<p 

The flat section of the potential does not guarantee that the solution of the field equations 
is of the form (1.2) and (1.3). The fact that the scalar (f>(t) rolls slowly corresponds to 
neglecting its first derivative (its "speed") in the Klein-Gordon equation (1.4), which 
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reduces to ~ — dV/ d<j). Alternatively, the slow-roll approximation corresponds to 

neglecting the kinetic energy density /2 in the expressions of the scalar field energy 
density and pressure 

P=y-V(0), (1-5) 

P = y + V(fl- (1-6) 

As a result, the scalar field is equivalent to a fluid with equation of state P ~ —p. For 
comparison, in de Sitter space the cosmological constant A can be regarded as a matter 
fluid with energy density and pressure 

PA= 8^' Pa = "8^' (L7) 

and equation of state Pa = — p\. 

Even if the potential V(<j)) has a flat section, the scalar <fi could still roll fast (with non- 
negligible 0) over it - the slow-roll approximation is an assumption about the solution 
(a(t) , <f>(t)) of the dynamical equations, not on the form of V ((/>). What makes this 
approximation viable is the fact that, in general relativity with a minimally coupled 
scalar field (and also when the field is non- minimally coupled [5, 6, 7]), de Sitter space is 
an attractor for the orbits of the solutions in phase space [4]. The main purpose of the 
present paper is to establish whether de Sitter space is an attractor also in more general 
gravity theories. 

There has been increasing interest in cosmology in alternative theories of gravity, 
with several different motivations [8]. One such motivation arises in the quest for a 
quantum theory of gravity: it is widely believed that quantum corrections modify the 
Einstein-Hilbert gravitational Lagrangian by adding terms proportional to higher order 
curvature invariants [9]. These corrections to the classical Lagrangian are small at small 
curvatures but become dominant when R grows, e.g., approaching a singularity. 

From another point of view, theories of gravity generalizing Einstein's relativity 
have been studied for decades at the classical level [8]. The prototypical alternative 
theory, Brans-Dicke theory, was originally motivated by the need to explicitly incorporate 
Mach's principle in relativistic cosmology, and has later been generalized to the class of 
scalar-tensor theories in which a Brans-Dicke-like scalar describes the gravitational field 
together with the metric tensor, and coupling functions appear in scalar-tensor gravity 
instead of coupling constants. In versions of these theories motivated by high energy 
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physics the Brans-Dicke-like scalar <fi is allowed to self-interact through a potential V(4>). 
The gravitational sector of scalar-tensor theories is described by the action 



Scalar-tensor gravity has been studied in relation to inflation, resulting in various ex- 
tended [10] and hyperextended [11] inflationary scenarios. Added interest comes from 
the fact that a gravitational scalar field is an essential ingredient of modern high energy 
theories unifying gravity with the other fundamental interactions (in particular string 
theories) [8], from certain similarities between scalar-tensor and string theories [8], and 
from the fact that the low-energy limit of the bosonic string theory is a Brans-Dicke 
theory with parameter uj = — 1 [12]. 

In this paper non-linear gravity and scalar-tensor theories are considered simultane- 
ously as special cases of the generalized gravity theory described by the action 



We neglect matter contributions to the action because we want to study situations in 
which the scalar <fi dominates the dynamics of the universe, such as during inflation 
in the early universe or in a late era in which a quintessence scalar field has come to 
dominate. The action contains the Ricci scalar R but no other curvature invariant. In 
addition to the simplification that accompanies it, this choice is motivated by the fact 
that in the homogeneous and isotropic cosmologies that we consider all the quadratic 
invariants of the Riemann tensor can be expressed in terms of R. The action (1.9) 
includes as special cases Brans-Dicke theory, scalar-tensor theories, induced gravity, 
quadratic Lagrangians, the theory of a scalar field coupled non-minimally to the the 
Ricci scalar, general relativity with or without a minimally coupled scalar field and a 
cosmo logical constant, and the theory of phantom fields. 

In the literature, slow-roll inflation in the context of generalized gravity is often 
considered. It is interesting to determine whether de Sitter spaces are actually solutions 
of the theory and whether they are attractors for the orbits of the solutions. This 
issue is crucial for understanding inflation in generalized gravity because the slow-roll 
approximation is meaningless unless there is a de Sitter attractor in phase space. 

Another independent motivation comes from the recent discovery [13] that the present 
expansion of the universe is accelerated, which has led cosmologists to postulate the ex- 
istence of a new form of of energy called quintessence or dark energy with the exotic 
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equation of state P < — p/3. Indeed, there are claims of evidence for a very negative 
pressure P < —p, a fact that, if confirmed, has interesting implications for the future 
of the universe - it could lead to a Big Rip singularity in a finite future [14] (see Sec. 6 
for a discussion). An obvious candidate for dark energy is the cosmological constant 
associated with de Sitter space. However, the cosmological constant carries with it two 
embarassing problems: 1) the well-known cosmological constant problem [15] of why the 
value of pa predicted by quantum field theory is 120 orders of magnitude larger than the 
energy density of the universe; and 2) the cosmic coincidence problem of why the dark 
energy is beginning to dominate the cosmic dynamics right now when there are galax- 
ies and human observers to notice it. These problems are only solved by an enormous 
amount of fine-tuning. For these reasons, theoretical models of quintessence explore 
different avenues. Among the many models proposed, modifications of Einstein gravity 
including non-linear corrections to the Einstein-Hilbert action have been proposed, in 
both the Einstein-Hilbert [16]-[23] and the Palatini form of the variational principle [24]- 
[26]. Such models do not usually admit a Minkowski solution that would be useful to 
study the weak-field limit of the theory - a de Sitter space is used instead for this pur- 
pose. Moreover, quintessence models that do not end in a Big Rip often evolve to a de 
Sitter phase in the future. Thus, both classes of models - either invoking a scalar field as 
dark energy (in general relativity or in scalar-tensor gravity), or advocating non-linear 
corrections to gravity, exhibit aspects related to the existence of de Sitter solutions. 

In a different context, it is interesting to examine the stability of general relativity 
with respect to small deviations from Einstein's theory due to quantum corrections. This 
is the approach adopted, e.g., in Ref. [27]. 

The purpose of the present paper is to establish conditions under which de Sitter 
solutions exist in the generalized theory described by the action (1.9), and to study 
their stability with respect to inhomogeneous perturbations. Our main motivation is to 
establish a firm foundation for the slow-roll approximation to de Sitter-like inflation in 
these theories. 

The issues of existence and stability of de Sitter solutions have been addressed in 
the literature only for special cases of the general theory (1.9) and usually only for 
spatially homogeneous perturbations. This limitation is probably due to the fact that 
inhomogeneous perturbations are in general gauge-dependent and they must be analyzed 
in the context of a covariant and gauge-invariant formalism. The latter is substantially 
more complicated than the analysis of time-dependent homogeneous perturbations. In 
the present paper the covariant and gauge-invariant formalism of Bardeen-Ellis-Bruni- 
Hwang-Vishniac is employed to study stability. This formalism has been used before to 
analyze the stability of de Sitter solutions in the theory of a scalar field coupled non- 
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minimally to the curvature [6, 7], and the stability of Einstein space in general relativity 
[28]. Following the same line of reasoning, it is also interesting to consider the stability 
of Minkowski space solutions of the theory. 

An independent motivation for the study of de Sitter space arises from the idea that 
the universe could have originated in a de Sitter state, thus avoiding the initial big bang 
singularity and evolving into an inflationary phase. Variations of this idea include the 
possibility of a Minkowski [29] -[32] or an Einstein space [33, 28] as a possible initial 
state. We include Minkowski space in our analysis as a special case of de Sitter space. 

The plan of this paper is as follows. In Sec. 2 we summarize the field equations 
of the generalized theory and we derive the conditions for the existence of de Sitter 
solutions. Section 3 addresses the issue of stability with respect to inhomogeneous 
perturbations using a covariant and gauge-invariant approach. Section 4 discusses the 
existence and stability of Minkowskian solutions of generalized gravity, while Sec. 5 
contains a discussion and the conclusions. We use units in which the speed of light 
c = 1 and 8ttG = 1, where G is Newton's constant, the metric signature is —,+,+, +, 
and □ = g^VaVb denotes d'Alembert's operator. For ease of comparison with previous 
works, the other conventions follow Refs. [34, 35]. 

2 Generalized gravity, fixed points, and de Sitter 
solutions 

Variation of the action (1.9) leads to the field equations of generalized gravity [35] 




g ab (RF — f + 2V) 



+V a V b F-g ab nF] , 



(2.1) 



1 (du n . ,„ , df dV\ 



(2.2) 



where 




l-Kr 



(2.4) 



(2.3) 
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in comoving coordinates (t,r,6,ip), the field equations assume the form 



H = --L(^ + F-HF) + ^, (2.6) 

where H = a/a is the Hubble parameter, an overdot denotes differentiation with respect 
to the comoving time t and the Ricci curvature is R = 6 (h + 2iJ 2 + K/a 2 ^. Only two 
equations in the set (2.5)-(2.7) are independent. 

There is now substantial evidence that the universe has flat spatial sections [36] and 
therefore from now on we restrict ourselves to the spatially flat case K = 0. In this case 
one can choose H and as dynamical variables - this is not possible if K ^ 0, in which 
case one has to consider as dynamical variable the scale factor a appearing in the field 
equations through the terms ±K/a 2 instead of H. However when K — these terms 
disappear and a appears only in the combination H = a/a and in its time derivatives 

(in non-linear gravity the field equations are of fourth order and H and H appear in 
the field equations). The phase space picture of the dynamical system depends on the 
specific form of the functions f((f),R), u{4>) and V(<f>). However, for any choice of /, uj 
and V, the fixed points of the system (if they exist) are given by 

(H, H, H, H, 0, 0) = (H , 0, 0, 0, 0o, 0) , (2.8) 

where H and O are constants, i.e., they are de Sitter spaces with constant scalar field 
[37]. The conditions for the existence of de Sitter fixed point solutions are obtained by 
substituting eq. (2.8) in eqs. (2.5) and (2.7), which yields the two conditions 

6# 2 F -/ + 2\/ = (2.9) 

and 

/o'-2Vb , = 0, (2.10) 



where 



F ° = M ' (2 ' n) 
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fo = /(0o,-Ro) , 



(2.12) 



V = V (0 O ) , 



(2.13) 



fo' 



dV 
~dUj> 

df 



4>o 



(4>o,Ro) 



(2.14) 
(2.15) 



and Rq = 12Hq. There are two independent conditions (2.9) and (2.10) for the existence 
of de Sitter solutions because only two equations in the set (2.5)-(2.7) are independent. 

Let us consider a few examples of specific gravity theories. Eq. (2.9) generalizes the 
condition 

6# 2 F -/ + 2A = (2.16) 
found in Ref. [27] for the non-linear gravity theories given by the choice 



1 , / = f(R) , V = A = const. 



(2.17) 



(there is only one condition in this case because the scalar field is not a dynamical 
variable). Note that not all generalized gravity theories admit de Sitter solutions. For 
example, non-linear theories with f(R) = AR n , A =constant, n > 2 and = 0, V = 
do not satisfy eq. (2.9) [27]. 

In general relativity with a cosmological constant A > and without scalar eq. (2.9) 
produces the familiar de Sitter solutions 



[H () . <>«) ■- ■ \ -.0 

If a minimally coupled scalar is present, de Sitter space is achieved if 

V >0, H = ± ] j^-, and V£ = . 
In the theory of a non-minimally coupled scalar field corresponding to 



(2.18) 



f{<P,R)=R{l-i<P 2 ) , u = l, 



(2.19) 



(2.20) 
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where £ is a dimensionless coupling constant, eqs. (2.9) and (2.10) reduce to the condi- 
tions for the existence of de Sitter fixed points previously found in Refs. [5]- [7] 



Hi (l - Ml) = | , (2.21) 

12£# 2 0o + VI = . (2.22) 
Recently the higher derivative theory of gravity described by 

ftf,R) = R- 0=1, u = 0, V = 0, (2.23) 

where /i _1//2 is a length scale, has attracted attention as a model theory for the accel- 
eration of the universe that does not require dark energy [16]- [26]. In this theory the 
correction to the Einstein-Hilbert Lagrangian is small at large curvatures, but becomes 
important as the universe expands and R — > 0. The conditions (2.9) and (2.10) for the 
existence of de Sitter solutions reduce to 




(2.24) 

Note that this theory does not admit a Minkowski space solution corresponding to H = 0. 



3 Perturbations 

If de Sitter fixed points exist for the dynamical system (2.5)-(2.7) with K — 0, the 
problem arises whether these fixed points are attractors in phase space or are unstable - 
a stability analysis is required to answer this question. The approaches to this problem 
available in the literature [38] -[40] [27, 6, 7, 41, 28] are limited to special cases of general- 
ized gravity theories and, usually, to homogeneous perturbations. The consideration of 
more general inhomogeneous perturbations is complicated by the gauge-dependence of 
this kind of cosmological perturbations. A gauge-independent analysis requires the use 
of a covariant and gauge-invariant formalism, which has been used before to study the 
stability of de Sitter solutions against inhomogeneous perturbations in the special case 
of the theory described by eq. (2.20) [6, 7]. Another problem addressed in the literature 
with a gauge-independent approach is the stability of the Einstein universe in general 
relativity with a non-minimally coupled scalar field [28]. 
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We proceed by using the covariant and gauge-invariant formalism of Bardeen [42, 3] 
further developed by Ellis, Bruni, Hwang and Vishniac [43, 44]. A version for generalized 
theories of gravity is given in Refs. [35, 34]. The metric perturbations are defined by 



#00 



1 + 2AY) , 



(3.1) 



9oi 



-a 2 BY, , 



(3.2) 



9ij = a 2 [hij (1 + 2H L ) + 2H T Yij] , (3.3) 
where the scalar harmonics Y are the eigenfunctions of the eigenvalue problem 

v i v i y = -fc 2 y. (3.4) 

Here is the three-dimensional metric of the FLRW background and the operator V* 
is the covariant derivative associated with h^, while k is an eigenvalue. The vector and 
tensor harmonics Y and Yy are defined by 



y 



Yi = ~ViY, 



(3.5) 
(3.6) 



We use Bardeen's [42] gauge-invariant potentials and $^4 and the Ellis-Bruni [43] 
variable A$ defined by 



H 



Ht a / a ■ 

H ^- + k{ B -k H - 



<$>A = A 



B 



En 



B 



1 



(3.7) 
(3.8) 

*~v (3 ' 9) 

Equations analogous to eq. (3.9) define the gauge-independent variables AF, Af, and 
AR. To first order the perturbations evolve according to the equations [35] 



A0 = 5(j)+r<j>(B -tH t 



y uo a<f> I 



k 2 



+ 



1 duo 



<j>($ A -3Q H ) + 



®a fdf 



CO 



, - , - 

2 d(p \ud<j> J 

2 dV\ 1 d 2 f 

d(f) J 2uj d(j)dR 



d_ (J_d£ _ ldV' 
d<f> y 2uj d(f) uj d(f) / 

AR , 



A<p 
(3.10) 
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AF + 3HAF+[ — --)AF + -AR+-u>(f)A<j)+ 11 



" ■ '2—^- ■■ \'^-) A0 



F ($ A - 3<5> H ) +-(FR-2f + AV) $ A , 



3\ d(f) d<fi d(j> ) 

(3.11) 



/ F\ ■ k 2 
H T + [3H + - )H T + — H T = , 



/ F \ 1 ( AF AF u ■ \ 

* + - *^ = - — -if— + -*a* , 



(3.12) 
(3.13) 



H 



lfu;- 2 3F 2 \ 1 f 3 FAF ( ■ k 2 3H F\ AF 

- I -0 2 + -— ) $a = - { -^r- + ( 3// 



2F 2 



2 2 F 2 



*#-4 + 2 # +6 " \ H+ 2F, 



AF 



a 2 2 FJ F 
A^j , (3.14) 

(3.15) 



<!>„ - //•!• f j H + — ) (2$ H - 4 — ( f 21 /// <!>,, 



F 

2F/ 



1 

2F 



AF orr AF ,„ AF u : A : 1 /; 2 dw 9/ n dV\ . / 



where AF = d(AF)/dt, etc., 
AR = 6 



+ 4# $ H + - if $ A - M2if + 4F, 2 - ^ J $ A 



and the effective energy density and pressure of the scalar are given by 



1 



— |— + ^ [RF — f + 2V) — 3HF + V C F C 



(3.16) 



(3.17) 



(3.18) 
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^ + l(f-RF-2V) + F + 2HF-^ V C F C 



(3.19) 



Here F c = h d c V d F is the spatial projection of the gradient of F. In the de Sitter 



background (2.8) the gauge-invariant variables reduce, to first order, to 

A0 = 50, AR = 5R , AF = 5F , Af = 5f, 
and the equations they obey reduce, to first order, to 

k 2 1 



A0 + 3# O A0 + 



a 2 2ujq 



A0 = ^ AF , 



AF + 3F AF + - 4F 2 J AF + y AF = , 



£2 

F^ + 3F -^r H — 5" Ht = , 



a; 

1 /AF 



2 V 
1 AF 



AF\ 



H 



2 F ' 
AF 



(3.20) 

(3.21) 

(3.22) 
(3.23) 
(3.24) 
(3.25) 
(3.26) 



$ H + 3H $ H - H $ A - 3H 2 $ A = -H ^ + ^l A/ 



2 F 



Fq 2 Fq 



with 

AF = 6 
and where 



2 k 2 ( k 2 \ 

§ H + 4F $h + - — $ H - F $4 + I — -4F 2 J $ A 



f <9 2 / 

J<t>R = 



d<pdR 



Jrr = 



d 2 f 



W>o,«o) 



OR 2 



,„ = d 2 f 

Jo — fli2 



(<t>o,Ro) 



(0o>-Ro) 



(3.27) 



(3.28) 



(3.29) 
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The comparison of eqs. (3.25) and (3.26) yields 



2F 

which, substituted in eq. (3.28), leads to 

' k 2 



AF 

$ ff = $ A = _ (3.30) 



AR = 6 



<$> H + ZH <$> H +\--ml)<$> H 



a 2 



(3.31) 



Eqs. (3.21) and (3.23) do not change form, while the remaining equations reduce to 
identities. The decoupling of scalar, vector and tensor modes is not apparent in the for- 
malism used, with the exception of tensor modes described by the perturbation H T . This 
quantity obeys eq. (3.23), which is decoupled from the other modes. We do not consider 
vector modes described by the quantity B since, as proven in Ref. [35], vorticity modes 
cannot be generated in generalized gravity when matter contributions are absent (i.e., 
when the scalar field or nonlinear corrections to R dominate). The vector mode B ef- 
fectively disappears from the gauge-invariant variables Acf> and §a = &h = — AF/(2F ) 
defined by mixing scalar modes and B. More naively, other authors refer to these facts 
by saying that the vector perturbations can be gauged away. For this reason in the 
following we consider explicitly only scalar and tensor perturbations. 

3.1 Stability with respect to tensor perturbations 

The evolution of tensor perturbations is regulated by eq. (3.23), where a(t) = a e Hot . 
By introducing the auxiliary variable 

u = aH T (3.32) 
and using conformal time rj defined by dt = a drj and the standard relation 

e Hot = 1— (3.33) 

a H r] 

valid in de Sitter space, eq. (3.23) is reduced to the formal Schrodinger equation 

S + (k 2 - 1) u( V ) = • (3.34) 



dr) 
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Let us consider the expanding (Hq > 0) de Sitter spaces (2.8). We are interested in 
the late time evolution of perturbations, corresponding to t — > +00 and r\ — > 0~. In this 
regime, the general solution of the asymptotic equation 



d 2 u 2 
drf r] 



, o«fa) = (3.35) 



is 



U (r,) = ?± + C 2 ri 2 (3.36) 
77 

for rj ^ 0, where C lj2 are integration constants. Then 

// T = -// (Ci + C 2 ?7 3 ) , (3.37) 

and the gauge-invariant tensor perturbation doesn't grow when 77 — > 0~. Hence, expand- 
ing de Sitter spaces are always stable with respect to tensor perturbations. 

Let us consider also the contracting (H < 0) de Sitter spaces (2.8). In this case 
t — > +00 corresponds to rj — > +00 and the asymptotic solutions of eq. (3.34) are free 
waves e ±tkri , hence the amplitude of the tensor perturbations 

±ikr] 

H T = - e |Ho|t (3.38) 

a 

diverges when t — > +00. As a conclusion, contracting de Sitter spaces are always unstable 
with respect to tensor perturbations. 



3.2 Stability with respect to scalar perturbations 

At a first glance it might seem that we are left with only one equation (3.21) to determine 
the perturbations A<fi and <3># = &a = — AF/(2F ), but this is not the case since one 
can Taylor-expand the coupling function / (0, R) obtaining 

^ = ^A0+^Ai?. (3.39) 

Eqs. (3.30) and (3.39) then yield 

AR=^$ H - |^A0, (3.40) 

J RR J RR 
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while eq. (3.21) becomes 



A0 + 3# O A0 + 



k 2 



2uq 



RR , 



A</> + 



H 



RR 



(3.41) 



The comparison of eqs. (3.40) and (3.28) yields 



$ H + 3H $ H +(%- Mil + ^- ) $h 



6/ 



A0 = 



RR 



(3.42) 



In the rest of this section we consider the case in which f RR ^ 0. This restriction leaves 
out linear theories of gravity, including general relativity, which are discussed in the next 
section. 

We have now the system (3.41) and (3.42) for A0 and which can be simplified 
by switching to the variables v and w defined by 



<5> H = 



w 



(3.43) 



and by using conformal time r\ instead of t. In terms of these new variables it is 



where 



d 2 v 



a 



P 



d 2 w 
drf 
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+ [k 2 + ^-\w + — v = 0, 



a = 



F 



P 

7 
S 



^IrrHq 

f<t>R 
SfR R H$ ' 

1 



6, 



2u 
F 



f2 ~ 

ot/" f" i J <t>R 
zv o Jo 



f. 



RR, 



^qHq Irr 



- 2 



(3.44) 
(3.45) 

(3.46) 
(3.47) 
(3.48) 
(3.49) 
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The system (3.44) and (3.45) can be linearized around the conformal time 7]q at which 
the perturbations originate and then rewritten as the first order system 



x 



w' = y , 



x 



y 



a 



-[k 2 + -z)v- , 



Vo 



— \ k 2 + ^7 }w 



This can be written in compact form as 

/ v' \ 



w 
x' 



M 



\ y' ) 



( V \ 

w 

X 

V y J 



where the matrix M is 



M = 



I 










i 


\ 















1 


- 


'r + 1 ) 


-1 
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i) 








(3.50) 
(3.51) 

(3.52) 

(3.53) 



(3.54) 



(3.55) 



A stability analysis for this system (see Appendix A) yields the result that de Sitter 
space is stable if any one of the following conditions is satisfied: 

• if bi > and C\ > 0, it must be k > k\ for stability 

• if b\ > and c\ < 0, it must be k > max {k±, k^} 

• if bi — and c\ > 0, it must be k > k\ 
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• if b\ = and c\ < 0, it must be k > max {k±, k±} 

• if -2 < bi < and c\ > bj/A, it must be k > h 

• if — 2 < 6i < and < c\ < 6f/4, it must be fci < /c < /c 5 (this inequality can be 
satisfied only if /c 5 > k\) or k > max{/ci, A; 6 } 

• if — 2 < < and C\ < 0, it must be k > max {ki, k 7 } 

• if bi < —2 and c x < 0, it must be k > max {k±, k 2 , k 7 } 

• if £»i < —2 and < c\ < bj/A, it must be max {k±, k 2 } < k < k 7 (when k 7 > k±, k 2 ) 
or k > max {k±, k 2 , k & } 

• if bi < —2 and c\ > bf/4, it must be k > max{fci, k 2 }. 

The parameters b\,c\ and the critical wave vectors ki are defined in terms of the 
values of the coupling functions and parameters by 



bi 



1 

7o 



1 ,« „// f^R 



2tUn 



2^ " - fS + 



+ 



Irr J 3/rrH 



- 2 



(3.56) 



Cl = 



2u 



-i 2 



8 + 



RR, 



3/rrH 



F 



1 



3JrrHq 2a;, 



+ 2V "-/5' + 



/, 2 



0.R 



+ 



<f>R 



(3.57) 



2u 



P 



+ 



FoP 



<f>R 



RR 



fe = lag 



2^n 



0R 



RR, 



3/rrHq 
F 



1/2 



3-f^o/i?R 



1/2 



(3.58) 



(3.59) 



fca = 



-h + J6f - 4ci 



(3.60) 
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c l 



1 1/4 



(3.61) 



k, 



5,6 



IftiliVftf-^x 



A; 7 = 



1 



|6i 



63 



Cl 



(3.62) 



(3.63) 



The effective Jeans wavelengths A; = 2tc/ fcj are determined by the values of f(<j>, R), V, 
and ui at the point (<f>o, Ro) and by the value of the initial conformal time 770 ■ The con- 
ditions for stability, when expressed in terms of f,V and lvq,t]o are n °t particularly 
illuminating - they become more meaningful when specific forms of these functions are 
adopted. Although the classification of the stability regions seems involved, once the 
generalized gravity theory is specified, the values of the parameters bi and c\ are com- 
pletely fixed and only one of the cases contemplated in the list above applies. Therefore 
this list contains complete conditions to check at a glance whether de Sitter space is 
stable against inhomogeneous perturbations in general theories of gravity. 

As an example consider the theory described by 



f = R + eR 2 



V = A > , 



(3.64) 



where e is a small positive constant (general relativity is recovered in the limit e — > 0). 
One has 

f = 12H 2 (l + 12e# 2 ) , F = 1 + 24e# 2 , V = A , (3.65) 
Vq =f' = ft = Ur = Sh = > f*R = 2e ( 3 - 66 ) 



and Ho = yA/3 as in general relativity. Since a = (2eA) 1 — 2 and (5 — 7 
eqs. (3.44) and (3.45) reduce, in the limit e — > 0, to 



v" + 



k 2 + 



1 



0, 



0, 



(3.67) 



w" + k 2 w = , (3.68) 

which constitute decoupled equations for the variables v and w with positive (angular) 
frequency squared and hence describe oscillatory perturbations associated with stability 
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of de Sitter space. The small correction eR 2 to the Einstein-Hilbert Lagrangian does 
not destabilize de Sitter space. Naively, this correction "reinforces" the effects of R. If, 
however, e were allowed to be negative, the effect of eR 2 would be the opposite and it 
would tend to destabilize de Sitter space. This is evident in eq. (3.67) when e < and 
the term containing e dominates for e — > 0, making the effective frequency squared of v 
negative and introducing exponential solutions associated with instability. 



4 Scalar-tensor theories 



We now restrict the stability analysis to scalar-tensor theories, for which the coupling 
/ ((/>, R) is linear in R and f RR = - this case was not contemplated in Sec. 3. Having 
already established that contracting de Sitter spaces are always unstable with respect to 
tensor perturbations, we only need to consider an expanding de Sitter background (2.8) 
with Hq > 0. It is straightforward to derive the first order evolution equations for the 
perturbations. Eqs. (3.30), (3.39) and (3.40) yield 



= - - 



1 f<t>R 



2 Fn 



A0 



while eqs. (3.31) and (4.1) lead to 



AR = - 



F 



A0 + 3H A0 




(4.1) 



(4.2) 



In conjunction with eq. (3.21) this yields 



A0 + 3H A0 + 



r,2 



fit 

1q y" 4- 

2 V + 



6/, 



— H 



UJQ 



x 2uoF J 



A0 = 



(4.3) 



if 1 + 3fl R /(2u F ) ^ 0. In the case in which 1 + 3f} R /(2w F ) = 0, instead, eq. (3.21) 
yields either the trivial solution A0 = or f ' — 2V Q " = 8lu Hi [46] . 

We look for solutions of the asymptotic form of eq. (4.3) at late times satisfying the 
ansatz 

A(j) = ee st , (4.4) 
where e and s are constants, with s satisfying the algebraic equation 



s 2 + 3H Q s + c = , 



(4.5) 
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with 



-3H ± A/9iJ 2 - 4c 
*± = 2 — |— 2 ( 4 - 7 ) 

of eq. (4.6) are such that i?e (s_) < 0, while the sign of Re (s+) depends on the sign of c. 
If c > then Re (s + ) < and there is stability. If instead c < then i?e (s + ) > and 
the de Sitter space (2.8) is unstable. The condition for the stability of de Sitter space 
in a scalar-tensor theory described by the action (1.8) is then 



c =-- 7 S f \ ■ ( 46 ) 



The roots 



f'd _ \rn _|_ 6 ^fl zt2 
2 K Fn -"0 



- i 1 + Si 



< . (4.8) 



• Non-minimally coupled scalar field 

As a particular case of scalar-tensor gravity we consider the theory of a non-minimally 
coupled scalar field given by the choice (2.20) of the coupling functions, yielding 

V IV 

^0 = 1-^8, /o = -2^o0o = 2V{ , U R = -2^ = , f» = -° , (4.9) 
where eq. (2.10) has been used. If the effective gravitational coupling of the theory 

is positive (which happens for any negative value of £ or, if £ > 0, for \(f>\ < (87rG£) ), 
then the denominator on the left hand side of eq. (4.8) is also positive and the stability 
condition of de Sitter space reduces to 

2 F 
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Upon use of eq. (4.9), this condition is written as 

Vo > /(*) ^ , (4.12) 

where 

x = £<fi, f(x) = T -^. (4.13) 

Eq. (4.12) coincides with the stability condition found in Refs. [6, 7] when O 7^ [45]. If 
the effective coupling (4.10) is instead negative, the stability condition is given by (4.12) 
with the direction reversed. 

In the case O = one has / = Rq = 12 H$, F = 1, f = J^r = 0, /„ = -2£R , and 
the condition for the stability of de Sitter space becomes 

U " + 12£# 2 > . (4.14) 

Using eq. (2.9), this assumes the form 

U " + 4£Uo>0, (4.15) 

the stability condition found in Ref. [6, 7] for cf) = 0. Eq. (4.8) generalizes to arbitrary 
scalar-tensor theories the stability conditions already known for non-minimally coupled 
scalar field theory, which are recovered as a special case. 



• General relativity 

In Einstein gravity with a minimally coupled scalar, 

u = F = l, fo = l2H 2 , ^ = /^ = /^ = / JiR = 0, (4.16) 

and the de Sitter space obtained if H = yVo/3, Vq = is stable if Vq > 0, in particular 
if the potential has a minimum at 0o- Hence in this case the concavity of the scalar field 
potential is the stabilizing factor, while its convexity would instead cause instability. 

If the scalar is absent the de Sitter space obtained thanks to a positive cosmological 
constant is automatically guaranteed to be stable with respect to inhomogeneous pertur- 
bations. It is well known that this space is also stable with respect to large anisotropic 
perturbations, with the exception of highly positively curved Bianchi IX models, as 
described by the cosmic no-hair theorems [47]. 
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• Phantom field 

The superstring-inspired theory of a phantom field with negative kinetic energy cor- 
responds to f((p,R) = R and u — — 1. The conditions for the existence of de Sitter 
solutions are 

Hi = | , V* = 0, (4.17) 

while the condition for stability reduces to V " < 0. Thus, de Sitter fixed points (H , O ) 
are attractors if Vq has a maximum at <f>Q. Due to the negative sign of its kinetic energy 
the phantom field "falls up" and settles in the maximum of the potential, de Sitter 
attractors have been found in super accelerating models of dark energy, thus avoiding 
evolution of the universe in a Big Rip singularity in the future [49, 50, 51]. 



5 Minkowski space and its stability 

Generalized gravity often admits Minkowski space solutions. In general, flat space so- 
lutions in generalized gravity are physically non-trivial and correspond to a balance 
between gravity and the scalar field <j) formally acting as a material source. In non- 
linear gravity theories where this balance cannot be achieved, a Minkowski solution may 
not exist. This is the case, for example, of the theory described by eq. (2.23). 

It has been suggested that the present universe could have originated from Minkowski 
space [30, 31] or from a static Einstein space [33, 28]. To pursue this idea it is necessary to 
ascertain the stability of Minkowski or Einstein space. The stability of Einstein spaces is 
a long-standing issue [38, 39, 40] and it has recently been revisited in general relativity 
by considering inhomogeneous and anisotropic perturbations [28]. Here we consider 
the stability of Minkowski space in generalized gravity with respect to inhomogeneous 
perturbations. 

The conditions for the existence of Minkowski solutions (with H = 0) of the field 
equations of generalized gravity are 

/o-2F = 0, (5.1) 

/o - 2V ' = . (5.2) 

A positive value of / , which describes a realistic situation, can be balanced in eq. (5.1) by 
a negative cosmological constant A, which is familiar to high energy physicists working 
with anti-de Sitter space. 
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In order to study the stability of Minkowski space one needs the linearized equations 
for the gauge-invariant variables 



A0 + 



A: 2 



1 



2uj 



A0 = I** AF , 
2lu 



AF + k AF + — AF = , 



Ht + k Ht = , 



(5.3) 

(5.4) 
(5.5) 



®H = -- 



1 AF 



2 F ' 
AF 



1 o 



Eq. (5.6) yields again $ H — -AF/(2F ) = $ A , while eq. (3.28) reduces to 

AF = 6 + k 2 . 



(5.6) 
(5.7) 

(5.8) 



Again, we consider scalar and tensor modes and drop the vector modes which cannot 
be generated in the absence of matter [35]. Inspection of eq. (5.5) allows one to con- 
clude at once that Minkowski space is always stable against tensor perturbations, which 
decouple from the other modes. To assess stability with respect to scalar perturbations, 
one considers the analogue of eqs. (3.41) and (3.42), which are 



^ 2 + r 2 + 3/l 



f<f>R n 
V + — W — , 



6/, 



RR 



d 2 W 



+ 



k 2 + 



2uj 



W - fs + p4>R 



RR, 



, Fo f<t>R n 
w H j— v — 

^0 JRR 



(5.9) 



(5.10) 



in the case fjm ^ 0, where the variables v and w introduced in eq. (3.43) now coincide 
with § H = §a and A0, respectively. 
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The system (5.9) and (5.10) can be reformulated as the first order system 



where 



In compact form, 



where 



x 



w = y , 



x 



- k 2 + 



3/ 



RR, 



6/ 



w , 



RR 



<*>o/ 



RR 



v + Dw , 
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k 2 + 



2u 



W - /o 7 + ^ 



/ 
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\y' J 
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\ y ) 
o 



V 



o 

f<t>B. 



1 0\ 
1 




^0 fRR 







(5.11) 
(5.12) 

(5.13) 
(5.14) 

(5.15) 
(5.16) 



(5.17) 



A stability analysis for this linear system, presented in Appendix B, leads to the result 
that Minkowski space is stable if 4c2 < b 2 , and one of the following conditions is satisfied: 

• b 2 = 0, c 2 < 0, and k > max {kg, k w } 

• b 2 < 0, C2 < 0, and k > max {kg, kis} 

• b 2 > and c 2 > 
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• b 2 > 0, c 2 < 0, and k > kg, 
where 



^ [ V ° h + fun) + 



c 2 = 



A: 8 



2uj, 

Fof^R 



3/ 



2 ' 



&9 

fcio 



A: 



13 



\ 



-b 2 


+ ^b\ - 4c 2 
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c 2 | 1/4 
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+ ^Jb 2 2 + A\c 2 \ 



(5.18) 
(5.19) 

(5.20) 

(5.21) 
(5.22) 

(5.23) 

(5.24) 
(5.25) 

Eqs. (3.39) and (3.30) yield AF = and <5> H = $4 = 0. Eq. (3.31) then yields AR = 
and eq. (5.3) becomes 

A0 + A; 2 A0 = O. (5.26) 

Eqs. (5.26) and (5.5) describe the evolution of scalar and tensor inhomogeneous pertur- 
bations of Minkowski space, which are effectively decoupled. It is obvious from these 
equations that Minkowski space is stable since the frequency squared is positive in each 
of these equations. 



Consider as examples general relativity and the theory described by 



f = R + eR 2 , 



1 , u = l, V = . 



General relativity is recovered by letting e — > 0. One has, for both theories, 
fo — , F — 1 , V — Vq — /q = /q = /^.r = = f RR = . 
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Let us consider now the class of scalar-tensor gravity theories, which have / (0, R) 
f((f))R and fun = 0. One has to consider eqs. (4.2) and (4.3), which become 

/o 

fo (2 V - f» 



AR = -^f(A^f> + k 2 A ( f ) ) , 



(5.27) 



A0 + 



k 2 + 



A0 = . 



(5.28) 



2/o + 3/^ 

The stability of a scalar perturbation depends on its wavelength and is achieved if 

1/2 

k>k u = J T, U , ;; ; u/ , (5.29) 



2/o + 3^ 

when the argument of the square root is positive (or for any wavelength if the latter is 
negative or zero). 



6 Discussion and conclusions 

Motivated by inflation, quintessence, and quantum gravity corrections to the low-energy 
gravitational action, we have derived conditions for the existence and linear stability 
of de Sitter solutions in a very general theory, including scalar-tensor gravity, induced 
gravity, non-linear gravity, 1/R corrections to the Einstein-Hilbert action, non-minimally 
coupled scalar field theory, general relativity with or without a minimally coupled scalar 
and a cosmo logical constant, and phantom fields. Our analysis does not depend on the 
specific form of the coupling functions and the conditions for the existence and stability 
of de Sitter solutions are, in this respect, very general. Minkowski space is studied as 
a special case of de Sitter space. The phase space picture of the theory depends in an 
essential way on the form of the coupling functions and the scalar field potential V((f>) 
but, in a spatially flat FLRW universe, the dynamical variables are always H and <fi. 
Although the field equations are of fourth order in non-linear theories of gravity, and 
hence the dimensionality of the phase space depends crucially on the form of f(<p,R), 
the fixed points of the dynamical system are always de Sitter spaces with constant scalar 
field. It is for this reason that the conditions for the existence of de Sitter space and for 
its stability can be expressed by inequalities valid for any choice of the coupling functions 
and the values of the free parameters. 

Eqs. (2.9) and (2.10) are necessary and sufficient conditions for the existence of de 
Sitter fixed points. Note that the existence of these de Sitter solutions is not auto- 
matically guaranteed in generalized gravity. Eqs. (2.9) and (2.10) reduce to conditions 
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previously obtained in special cases of generalized gravity (non-minimally coupled scalar 
field theory [6, 7] or theories with / = AR n [27]). 

When de Sitter fixed points exist, their stability against linear inhomogeneous per- 
turbations and their attractor behaviour are assessed by using a covariant and gauge- 
invariant formalism originally developed to study perturbations of FLRW spaces. It is 
established that expanding (resp. contracting) de Sitter spaces are always stable (resp. 
unstable) with respect to tensor perturbations. Scalar perturbations may threaten the 
stability of expanding de Sitter spaces, which are the ones of interest for inflationary 
and quintessence scenarios of the real universe. For non-linear theories of gravity (with 
fun 7^ 0) Subsection 3.2 provides the desired stability conditions, while section 5 de- 
scribes the stability of Minkowski space. For linear gravity theories with Jrr = 0, 
including scalar-tensor gravity and general relativity, the stability conditions for de Sit- 
ter space are given by eq. (4.8), while the stability of Minkowski space is determined 
by eq. (5.28). In the particular case of non-minimally coupled scalar field theory, these 
conditions reproduce those already known from a previous analysis [6, 7]. 

The analysis presented here can be generalized further. First, it is well known that 
there can be attractors in phase space that are inflationary but are not fixed points. 
This is the case of power-law inflation a(t) = a t p with p > 1, which is an attractor 
solution of Brans-Dicke cosmology when the only form of matter present is a cosmological 
constant. Extended inflationary scenarios [10] are based on the presence of this attractor. 
Power-law inflation is also an attractor in scalar-tensor theories generalizing Brans-Dicke 
gravity [8] and in many non-linear theories [48]. Hence, in general, the fixed points do 
not provide the complete phase space picture. Second, we restricted our attention to 
inhomogeneous perturbations. Although more general than the case of homogeneous 
perturbations usually studied in the literature, it would be interesting to generalize the 
stability analysis to anisotropic and to nonlinear perturbations. Finally, we considered 
only four spacetime dimensions but quantum gravity extensions of general relativity 
would call for a more general analysis in arbitrary spacetime dimension. 

The stability conditions derived here can be applied to the investigation of the Big 
Rip singularity in the future. It has been pointed out that the present expansion of 
the universe may be superaccelerated, i.e., H > 0, which is equivalent to an effective 
equation of state parameter w = P/p < — 1 for the dark energy dominating the dynamics 
of the universe at redshifts z < 1 [14]. Superacceleration cannot be achieved with 
a canonical, minimally coupled scalar field in Einstein gravity [8, 14]. If the universe 
really superaccelerates (which is not yet established due to the error in the observational 
determination of the parameter w) it runs the risk of ending in a Big Rip singularity 
in a finite future [14]. This kind of singularity is different from the Big Bang or the 
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Big Crunch because the universe expands explosively while the energy density of dark 
energy diverges instead of getting diluted, due to the peculiar equation of state P < —p 
[8, 14]. If the equation of state of the dark energy is constant with w =const.< — 1, 
the Big Rip is unavoidable. However, a time-dependent effective equation of state with 
w = w(t) is more realistic and in this case scenarios have been proposed in which the Big 
Rip is avoided. At present, there are in the literature superaccelerating models in which 
the Big Rip is unavoidable [14] and others in which a late time de Sitter attractor with 
H — exists which stops superacceleration and avoids the Big Rip [49, 50, 51]. It is in 
this context that the stability conditions derived here can play a role: these conditions 
help assessing the stability of de Sitter spaces and deciding whether a late time de Sitter 
attractor exists that attracts the orbits of the solutions of the field equations in phase 
space, thus avoiding the Big Rip. A generic statement about the fate of the universe 
requires the knowledge of the attraction basin of an attractor, and this issue can only be 
addressed in a specific theory with the form of the potential V(<f>) fixed. However, the 
conditions for the existence of de Sitter attractors provide an answer about the possibility 
of avoiding the Big Rip at least for initial conditions lying in a certain attraction basin 
to be determined. 

The stability conditions derived in this paper will be applied elsewhere to specific 
models of inflation and dark energy. 
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Appendix A: stability analysis for de Sitter space 



The stability of the system (3.54) is assessed by studying the sign of the real part of the 
eigenvalues A of the matrix M. The characteristic equation Det (m — \Pj , where / is 
the identity matrix, reduces to 



where 



A 4 + BiA 2 + Ci = 0, (A.l) 



B l = 2k' 2 + ^ + ^, (A.2) 



2k 2 + 


a 




Vo 


(W 


a 




Vo 



c> = i • : ( i> :1 • 3 ) - ^ • ( A -3) 



Vo) Vo 



(A.4) 



The squares of the roots are given by 



A| = , A x = B'l - AC, . (A.5) 

Let us consider the cases C\ > 0, C\ — 0, and C\ < separately. 

• < d < Bl/A 

Note that one cannot have B\ — in this case. If B\ > 0, then A^_ < and the 
roots A±± of eq. (A.l) are purely imaginary, i.e., Re (A) = and de Sitter space is 
neutrally stable. 

If Bi < then A^_ > and all the four eigenvalues A±± are real, with two of them 
positive and two negative. The positive ones give rise to instability. 

• Ci > B 2 JA (Ai < 0) 
If B l ^ 0, then 



-B 1 ±iJ\A, 



and 



K = f =pe w± (A.6) 



A ±± = ±v / Pe 4 -. (A.7) 



28 



In this case two roots A have positive real part and de Sitter space is unstable. 
If Bi = and C\ > one has A± = ±iy/d and 



A 



1/4 



V2 



;i±0 . 



(A.8) 



The roots with positive real part are associated with instability of de Sitter space. 

• d = 

In this case eq. (A.l) gives A^ = or A| 3 = —B\. If i?i > the roots A2,3 are purely 
imaginary, corresponding to oscillating perturbations and to stability. If instead B\ < 
there is a real positive root associated with instability. If B ± — d — then A = and 
de Sitter space is neutrally stable. 

• d < 



In this case 



A 



-B x ±^Bl + A\d 



(A.9) 



The lower sign produces two imaginary roots, while the upper sign gives two real positive 
roots associated with instability. 

As a summary, de Sitter space is stable if 



< d < 



and 



Bi > 



(A.10) 



and unstable otherwise. 

The conditions for stability can be formulated in terms of effective Jeans wavelengths 
and of the values of the coupling functions and parameters of the theory. The inequality 
Ci > is equivalent to 

2k 2 ril + a + 7 -(36 > , (A.ll) 



which can be expressed as 



[2rjl 



1 1 ( fit r,irll f4>R , FoftR \ , o 

Vo l 2uJ o V Jrr 6H oJrr/ 



Fn 



1 ^ V2 



"SJrrHo 



(A.12) 
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The inequality B\ > is equivalent to 

2k 2 r]% + a + 7 > , 



or 



where 



b x + 2 + A; 2 > , 
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(A.13) 
(A. 14) 
(A.15) 



fun J S/rrHq 

Eq. (A. 14) is satisfied for any k if 61 + 2 > or, if &i < —2, by wave vectors k such that 
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The inequality Ci < 5 2 /4 is equivalent to 
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which can be written as 
where 
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(A.19) 



To identify the values of the wave vector k that satisfy the inequality (A.18) one studies 
the sign of (p(k) by distinguishing several cases. 

• 61 > and ci > 0: the inequality (A.18) is satisfied for any value of k. 

• b± > and c\ < 0: the curve representing (fi(k) starts negative at k = and is 
always increasing, crossing the /c-axis at a point k 3 . The function (p(k) becomes 
positive for k > £j 3 , where 



-61 + \ b{ - 4ci 



(A.20) 
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• bi — 0: then, if c\ > 0, (p(k) > for any value of k > 0; 
if C\ < then </?(&;) > for k > A; 4 = Icij 1 ^ 4 . 

• &i < 0: the curve representing (p(k) starts from the value c± at k — 0, decreases 
for < k < \J~\b~i\ /2, reaches db liiiiiiiiiuiii ctt k — \J~\b~i\ /2, and then is always 



increasing for fc > y |&i| /2. If ci > and the minimum is non-negative, then 
<^(fc) > for any value of k and the equation <f(k) = has no real roots. This 
happens if c\ > bf/4. 

If instead C\ > but the minimum of <p(k) is negative there are two real roots 
of the equation (p(k) = and ip is positive for < k < k 5 and for k > k 6 , and 
negative otherwise. This situation occurs if c\ < b\/A and the critical wave vectors 
are 



k 



5,6 



|6i|± J6?-4ci 

• (A.21) 

Finally, if C\ < 0, the curve representing (p(k) is negative for < k < k 7 and 
positive for 



|&i| + \ b{ + 4 ci 
k>k 7 = ^ 111 -. (A.22) 

By putting together the three conditions for stability C > 0, B > and C < 
B 2 /4 one obtains that, in order for de Sitter space to be stable with respect to scalar 
perturbations, the latter must have wave vectors in one of the following intervals: 

• if bi > and C\ > 0, it must be k > k\ for stability 

• if b\ > and c\ < 0, it must be k > max {fci, /C3} 

• if bi = and ci > 0, it must be k > k\ 

• if 61 = and c\ < 0, it must be k > max {ki, k±} 

• if -2 < bi < and ci > 6f/4, it must be k > h 

• if —2 < 61 < and < c\ < £>f/4, it must be k\ < k < k 5 (note that this inequality 
can be satisfied only if fc 5 > ki, ^2) or > max {/ci, /c 6 } 

• if — 2 < £»! < and c\ < 0, it must be k > max {k 1: k 7 } 
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• if bi < —2 and c\ < 0, it must be k > max {k±, k 2 , kr} 

• if bi < —2 and < C\ < bj/A, it must be max{/ci,/c 2 } < k < k 7 or k > 
max {ki, k 2 , k e } 

• if bi < —2 and c\ > bf/A, it must be k > max{ki, k 2 }- 

Appendix B: stability analysis for Minkowski space 

The stability of the system (5.16) is determined by the real part of the eigenvalues A of 
N. The characteristic equation Det(iV — A/) is 

A 4 + J B 2 A 2 + C 2 = 0, (B.l) 

where 



B 2 = k 2 + ^^-D, (B.2) 



2 Fof^n ( f^R 



3 fun V 2w fun j 



The squares of the roots are given by 

A 2 ± = - B *Y*~ 2 , A 2 = Bl - AC 2 . (B.4) 
Let us consider the cases C 2 > 0, C 2 = and C 2 < separately 
• < C 2 < Bl/A 

Note that one cannot have B 2 = in this case. If B 2 > 0, then A^_ < 0, the roots 
A±± of eq. (B.l) are purely imaginary, and Minkowski space is neutrally stable. 

If i?2 < then A^_ > and all the four eigenvalues are real, with two of them positive 
and two negative - the positive ones make Minkowski space unstable. 



C 2 > Bl/A 
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If B 2 ^ 0, then _ 

-B 2 ± iJ~A 2 

A| = / =P*"* (B.5) 

and fl 

A±± = ±v/pe^. (B.6) 

In this case two roots A±± have positive real part and Minkowski space is unstable. 

If B 2 = and C 2 > one has A±± = ±2~ 1/2 C 1/4 (1 ±i). Two roots have positive 
real part, corresponding to instability. 

• C 2 = 

In this case eq. (B.l) gives A 2 = or A 2 = —B 2 . If B 2 > the roots are purely 
imaginary, corresponding to oscillating perturbations and to stability. If instead B 2 < 
there is a real positive root associated with instability. If B 2 = C 2 = then A = and 
Minkowski space is neutrally stable. 

• C 2 < 

In this case 

^-Wprg (B7) 

The lower sign produces two imaginary roots, while the upper sign gives two real roots, 
one of which is positive and is associated with instability. 
To summarize, Minkowski space is stable if 

B 2 

< C 2 < -1 and B 2 > (B.8) 
4 

and unstable otherwise. 

Let us express the conditions above in terms of the coupling functions and of the 
wave vector k. The inequality C 2 > is equivalent to 

ip(k) = k 4 + b 2 k 2 + c 2 > , (B.9) 



where 



\ | 2y" _ j" _|_ Ur \ _|_ FpftpR ^ 
2^o \ /rr/ 3f RR 
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The inequality B 2 > is equivalent to 



W - f») . 



(B.11) 




>0, 



(B.12) 



which is always satisfied if 62 > 0, and is satisfied only for perturbations with wave 
vectors such that 



Note that in this case the wave vector drops out and this is a requirement on the theory 
of gravity independent of the wavelength of the inhomogeneous perturbation. Next, one 
studies the sign of ip(k) as in Appendix A, with the following result. 

• b 2 > and c 2 > 0: then ip(k) > for any value of k. 

• b 2 > and c 2 < 0: it is ip(k) > for 




(B.13) 



when b 2 < 0. The inequality C 2 < B\ 



/4 is equivalent to 

4c 2 < b\ . 



(B.14) 




(B.15) 



• b 2 < 

If < 



• b 2 = 0: 

if c 2 < 





(B.16) 



Finally, if c 2 < 0, ip{k) is positive for 




(B.17) 
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The stability of Minkowski space is assured by imposing that the three inequalities 
(B.9), (B.12) and (B.14) hold simultaneously. It must be 4c 2 < b%, plus one of the 
following conditions must hold: 

• b 2 = 0, c 2 < 0, and k > max{/c 8 , k w } 

• b 2 < 0, c 2 < 0, and k > max{fc 8 , A; 13 } 

• b 2 > and < c 2 < b 2 2 /4 

• b 2 > 0, c 2 < 0, and k > kg. 
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